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Abstract 



We set up a suitable renormalization programme for the one-loop computation of the decay 
rate T{H -^ hh) in the Minimal Supersymmetric extension of the Standard Model. We 
then perform an explicit diagrammatic calculation, including the full contributions from 
top, bottom, stop and sbottom loops. We find that, for tan/3 close to 1 and mu ^ 2mj, 
,^ . the results can significantly differ from those previously obtained in the eff'ective potential 

approach. However, the latter method remains a good approximation in the region of 
parameter space which is most relevant for H searches at large hadron colliders. 



CERN-TH.6603/92 
July 1992 



^Also at Dip. di Fisica, Univ. di Padova, and INFN, Sezione di Padova, Padua, Italy. 
•^On leave from INFN, Sezione di Padova, Padua, Italy. 



Radiative corrections to the parameters of the Higgs boson sector in the Minimal Supersym- 
metric extension of the Standard Model (MSSM) have recently received much attention. 

After the discovery [|l| that top and stop loops can cause large corrections to the masses of the 
neutral CP-even Higgs boson^, radiative corrections to Higgs boson masses have been computed 
by a variety of methods: the renormalization group approach [^ , the effective potential approach 
[^^ and the diagrammatic approach 0-0. 

The renormalization group approach assumes that there are two (or more) widely separated 
mass scales, for example 

MsusY (~ m-ti ~ "^ta ~ • • • ~ ^H ~ mH± ~ ^a) > mz (~ mh ~ mt) , (1) 

and considers the effective theory for the degrees of freedom lighter than Msusy- It then solves 

(non-supersymmetric) renormalization group equations to obtain running parameters down to 

the scale Q = mz, imposing the tree-level relations of the MSSM as boundary conditions at 

the scale Q = Msusy- This approach has the advantage of resumming the leading corrections, 

proportional to logi^MsusY/^TT-z), so that even the case of Msusy orders of magnitude larger than 

mz can be dealt with in perturbation theory. On the other hand, if supersymmetry is to solve 

the naturalness problem of the Standard Model, one expects the various mass parameters of 

1/9 

the MSSM to be scattered around the electroweak scale, Gp ~ 250 GeV, so that assumption 

(|T]) breaks down. 

The effective potential approach consists in identifying the Higgs boson masses and self- 
couplings with the corresponding derivatives of the one-loop effective potential, evaluated at 
the minimum. By definition, this approach evaluates all Higgs self-energies and vertices at van- 
ishing external momentum. In the case of radiative corrections to the Higgs boson masses, this 
was shown to be a rather accurate approximation P,^. Actually, when the external momentum 
(i.e. the Higgs mass) approaches or exceeds the threshold of the internal particles, the full 
correction can be rather different from the zero-momentum one. However, in that case correc- 
tions themselves are small, either in the absolute sense or relatively to the (increased) tree-level 
mass. Other possible drawbacks of the effective potential approach are the gauge- and scale- 
dependence of the associated quantities^. These are not serious problems in the computation 
of the mass corrections: the dominant ones come from quark and squark loops, which do not 
introduce spurious dependences on the gauge parameter into the results; also, wave-function 
renormalization effects, responsible for the scale dependence, are generally small with respect 
to the overall mass corrections. 

The diagrammatic approach consists in performing the complete one-loop renormalization 
programme, specifying unambiguously the input parameters and the relations between renor- 
malized parameters and physical quantities. This approach gives the most precise computational 
tool in the case of supersymmetric particle masses spread around the electroweak scale, and re- 
sults which are formally gauge- and scale- independent. Since corrections can be numerically 
large, however, one has to pay attention and improve conveniently the naive one-loop calcu- 
lations when necessary. An example is the determination of the neutral CP-even masses, as 
discussed in detail in ref. 01. 



^Previous studies [EJ either neglected the case of a heavy top quark, or concentrated on the violations of the 
neutral Higgs mass sum rule, without computing corrections to the individual Higgs masses. 

^We recall that, although the effective potential is scale- independent, scale dependence enters its derivatives 
through the renormalized fields. More generally, the issue involved here is the dependence on the renormalization 
scheme, scale dependence being interpretable as a particular kind of scheme dependence. 



Whilst radiative corrections to Higgs boson masses are by now well under control, the study 
of radiative corrections to Higgs boson couplings is still at a less refined stage. In most phe- 
nomenological P,^,p!0|- [l^] and experimental [|1^ studies, radiative corrections to the Higgs cou- 



plings to vector bosons and fermions have been taken into account only approximately, by 
improving the tree-level formulae with one-loop corrected values of the h-H mixing angle, a, 
and with running fermion masses, evaluated at the typical scale Q of the process under consid- 
eration. Residual corrections are expected to be numerically small, with the possible exception 
of important threshold effects |I^Jl5[] . 

In the case of the Higgs boson self-couplings, which control decays like H — > hh, H — > AA 
and h -^ AA when the latter are kinematically allowed, it is known |^ that radiative corrections 
can be numerically large. Radiative corrections to cubic Higgs boson self-couplings have been 
computed, at different levels of approximation, both in the effective potential approach |[5|,p!6|,p 



and in the renormalization group approach [|I^. Given the fact that, in addition to the masses 



of the virtual particles in the one-loop diagrams, two different mass scales are involved in the 
decays H -^ hh, H -^ AA and h -^ AA, the mass of the decaying particle and the mass of the 
decay products, one might suspect that momentum-dependent effects, which are neglected in 
the renormalization group and in the effective potential approaches, could play a role. 

The purpose of the present work is to perform the diagrammatic calculation of r{H —>■ hh) 
at the one-loop level, and to compare the results with those obtained in other approaches. The 
main motivations for choosing this particular decay are the relative simplicity of the calculation 
and the fact that, even after the inclusion of the leading radiative corrections, H ^ hh is 
the dominant H decay mode over a large region of parameter space. A detailed discussion of 
the MSSM Higgs branching ratios at the one-loop level would require extending the present 
calculation to other decay modes, as currently under study. 

This paper is organized as follows. We begin by setting up a convenient one-loop renormal- 
ization programme for the decay rate r{H -^ hh) in the MSSM. We then perform a complete 
computation of the contributions due to top, stop, bottom and sbottom loops. Finally, we com- 
pare our results with those obtained at the tree level and in the effective potential approach, 
and we discuss their phenomenological implications for the detection of H at future colliders. 

The notation of the present paper will closely follow that of P,^, unless otherwise stated. In 
[^, radiative corrections to the neutral Higgs boson masses were computed in the DR scheme 
[p!8| , using the physical mass niA and /3 = [3^^{Q = mz) as input parameters, and explicit formu- 



lae for the physical masses rrih and mu were given. Here we shall adopt the same renormalization 
scheme for the computation of T{H -^ hh) at the one- loop level. We define the H and h fields 
in the CP-even neutral Higgs sector by if = cos a 5*1 + sin a S2 and h = — sin a Si + cos a S2. 
As explained in [^], for a satisfactory convergence of the perturbative expansion it is important 
to define the mixing angle a in terms of a mass matrix which includes the leading, momentum- 
independent one-loop self-energy corrections. 

The decay rate for the process under consideration reads 



r(ff^ftft)=l-^'^ -"")!' 



where the one-loop-corrected amplitude A is given by 



^< (2) 



A-'^iZjj )Hi{Z^ )hj{Zf^ )hkKjk + ^Hhh ■ (3) 



In eq. (0), i,j, k = h,H and Xijk are the tree- level cubic Higgs couplings in the CP-even sector| 

^HHH = T— 3 cos(/^ + a) cos(2a) , (4) 

z cos uw 

XuHh = —T, 7 — [—2 cos(/5 + a) sin(2a) — sin(/3 + a) cos(2a)] , (5) 

2 cos Uy\/ 

Xhhh = 7, — [2 sin(/? + a) sin(2a) — cos(/? + a) cos(2q;)] , (6) 

2 cos Ow 

Xhhh = -7^ ^— 3 sin(/? + a) cos(2a) , (7) 

2 cos t>w 



where g, mz = y{g'^ + fi''^)('^i +'^2)/2 and cos6'ty = gj \Jg^ + g''^ are (-D-R) renormalized pa- 
rameters. The last term Kjihh is the [DK) renormalized one-loop proper vertex, evaluated with 
on-shell external momenta. 

Finally, the 2x2 matrices Zn and Zh correspond to a finite wave function renormalization, 
which must be taken into account in the computation of a physical amplitude. Such corrective 
factors can be interpreted as the effect of inserting 'bubbles' on the external legs of the tree- 
diagrams. More precisely, the two matrices Z^ and Zh are nothing else than the (matrix-) 
residues of the renormalized (matrix-) propagator Gij{p^) of the neutral CP-even Higgses, at 
the poles p'^ = mj^ and p"^ = ml, respectively. Actually, Zh and Zh are extracted from the 
real part of the propagator, so in the following we shall implicitly consider only real parts of 
propagators and self-energies. The renormalized inverse propagator r(p^) was computed, in the 
(5'i,S'2) basis, in ref. P], and the one-loop-corrected masses niH and nih were extracted from 
it. We now rewrite r(p^) in the {H, h) basis in terms of mn and rrih: 

''^^ ' [ ThHip') p'-ml + Uhhip') - Uhhiml) ) ■ ^'^ 



In eq. (|^), the I[ij{p ) are one- loop Di^-renormahzed self-energies, and 

THhip^) = ThHip^) = ^sm2{P + a)flzz{ml) + ^sm2(P - a)AflAA{ml) + AUnhip^) , (9) 

where for a generic self-energy Il{p^) we define An(p^) = n(p^) — 11(0). Around the pole 
p^ — rriff, we have G(p^) ~ Zu/ip^ — ^h)^ where 

7 -( l-nW("^lf) -^Hhiml)/iml-ml)\ , . 

^ [-T^h{ml)/{ml-ml) ' ^'"^^ 



Around the pole p'^ = ml, we have G{p'^) ~ Zhjip^ — fnl), where 

7 ^ / -THhimD/iml-ml) 

^' [ -THhimD/iml -ml) 1 - Ii;,,(m^) 

In eqs. (|TUp and (^, il\p^) = dll{p'^)/dp\ 
Expanding eq. (^, we finally obtain 



:il) 



A 



-flnHi^n) - Khi^h) 



'2VHh{,mj^)XHHh-'^Hh{mH)Xhhh , . ..^.n, 

^Hhh-\ 9 9 \--^^Hhh- [i-Z) 



mtr — TTlt 

■^In the Feynman rules these couphngs are multiphed by i. Also, we denote by lA the proper vertex obtained 
from the Feynman diagrams. 



We now proceed to the evaluation of the general formula for A, eq. ([T2D, in the particular 
case in which only diagrams corresponding to top-stop-bottom-sbottom loops are taken into 
account. This should give the dominant one-loop correction to r{H — >■ hh) fn\ . 



The expressions for Il'jjfj{m'jj) and 11'^ 



hh 



,mt 



can be trivially obtained from the expressions 
of AIIhh and Allhh given in ref. p|. To obtain an explicit expression for THhip"^), eq. (|^), we 



need n^^(m|), Aflyi^(m^) and AYlnhip"^)- The first two were givenQ in ref. 
one we find: 



For the last 



AflHhip' 



Sg^m'f 



sm a cos a 



IGTc'^m'yy 



sin P 



+ 3m'fAF{mt, rrit, p) - 3G(mj, nit, p) 



Sg'^ml sin a cos a 



IGTT^m^ 



w 



cos^/5 



— + 3mlAF{mb, mi,, p) - 3G{mi„ nih, p) 
o 



^g' 



^Q^^^W g=t,ba,b=l,2 



Y. Yl (^Hq,qbChq,qAFirnq^,mq^,p) . 



(13) 



In eq. (^3|), p = \fpp' and AF{m\,m2,p^ = F{mi,m2,p) — F{mi,m2,0), whilst the coefficients 
correspond to the trilinear Higgs-squark-squark couplings and are summarized in the 



^in. 



naQb 



Appendix. The functions F and G were given in [Q. 

Finally, we need to compute Anhh- The three basic topologies of the diagrams contributing 
to Afjhh are depicted in fig. 1. Accordingly, the result for A^hh can be written as the sum of 
three contributions: 



where 



A 



Hhh 



A 



in 

Hhh 



A 



in) 

Hhh 



A 



[III) 

Hhh y 



(14) 



A 



(I) 

Hhh 



A 



{II) 

Hhh 



A 



{III) 

Hhh 



"ig^ml 



sm a cos a 



IQlt'^lTLy^ 



siif (3 



F^irrit, mt, mu) + 2F''{mt, rrit, rrih) 



+ Uml - -mjj - mlj f{mt, rut, mt, m^, mh, mh) 



Sg^m^l cos a sin a 
IGyr^mu/ cos^ P 



+ 

+ Umt 
2-3g' 

-3g' 



F^mh, mi,, mn) + 2F''{mb, m^,, my,) 



,nl-,nl f{,n,.,n,,rn,,,n„,rn„m,) 



Yl Y (^HlaQt ChqaQc ChMc f i^Qa , "^96 , ^Qc ]mH,mh,mh) 



(15) 
(16) 



1«^2™3 Z^ Z^[(^HqaqbChhqaqtF 
iOTT m^r q^j ,, ^ j, 

+ 2chq,g^CHhq^q,F\mq^, m,,, mh)] 



m, 



Qa 5 '^I^ljb ' 



.mn) 



(17) 



The expressions for the functions F'^ and /, as well as for the coefficients Cijq^q^ corresponding 
to the quartic Higgs-Higgs-squark-squark couplings, are collected in the Appendix. 



^The expression of Ilzz{rn^), which was correctly given in the preprint version, contains a misprint in the 
pubhshed version: the coefhcient appearing in the third hne should read (sfctL — c^CtR)'^ and not (sfctL + c^CtR)'^- 
There are other misprints in the published version, which were not present in the preprint version: in the 
conventions stated in footnote number 2, the part of a gauge-boson self-energy diagram proportional to g^i, 
should read —ig^^Il{p^); in the expression for A'j^j^/vn^, d*2 should read 2d\2', the remaining ones are obvious. 



This concludes our analytical evaluation of the one-loop-corrected amplitude and decay 
width, eqs. (^ and (|^). As a check equivalent to divergence cancellation, we have explicitly 
verified that the amplitude, and consequently the width, do not depend on the renormalization 
scale Q, as expected for physical quantities. Actually, such Q independence holds up to higher- 
order terms, consistently with the one-loop accuracy of our computation. It is the result of 
a cancellation between the explicit Q dependence of self-energies and proper vertices, and the 
implicit Q dependence of the parameters contained in Xuhh'- in our case, consistency requires 
that we consider only the Q dependence of parameters that is due to top, stop, bottom and 
sbottom virtual effects. In the following numerical calculations, we set Q = mz- 

Consistency would also require that we specify the input parameters g, mz and cos 9w [or, 
equivalently, a = g^ sin^ 6'vi//(47r), mz and cos 9^] in the DR scheme, and in a theory containing 
the stop and sbottom degrees of freedom besides the Standard Model particles. The DR mass 
m,z is related to the physical mass m,z,phys by 

^l = m^phys + ^zz{ml) . (18) 



Similarly, the DR fine structure constant a is related to a^^ ~ 137 by 

- = — + 7r- L ^i log , 19) 

where Aaught is the (large) contribution from charged leptons and the five observed quarks, 
and the remaining (small) contributions from the heavy particles in the model, denoted by the 
index i, are proportional to the corresponding one- loop QED /5-function coefficient bi. Finally, 
the DR electroweak mixing angle is given by 



2/1 ^^W,phys 1 , ^Wwi^T^w) ^Zz{mz) 

COS Uv' ~ 



'^ ~ 772 



^^Z,phys 



1 + 



m^ m? 



(20) 



At the level of accuracy of our numerical examples, however, it is enough to work with the 
fixed input parameters m^ = 91 GeV, sin^ 9w = 0.23 and a~^ = 128. This approximation 
is justified by the fact that there are other effects not accounted for in our results: one-loop 
corrections involving loops of gauge bosons, Higgs bosons, gauginos and higgsinos; also, order 
h"^ or g'^ (two-loop) corrections to the one-loop diagrams considered here. 

Before moving to the numerical evaluation of our results, we would like to relate them to 
the results one obtains in the effective potential approach, which consists in approximating the 
amplitude A by 

.e.p. ^ _ ( d'Veff \ 

\ / mm 

Since this amounts to computing the Hhh 3-point function at vanishing external momenta, we 
can easily obtain an explicit expression for A^''^' by taking a special limit of our previous result: 

A^'-'P- = Xuhh + J^Hhh\mt,=0,m„=0 ' (22) 



As a check, we have computed A^'^' from its definition (|2T|) and verified that eq. (0) gives an 
identical result. 

To illustrate our results, we show in fig. 2 the one-loop-corrected width T{H —>■ hh), as a 
function of mn, corresponding to four representative parameter choices. For simplicity, in our 



numerical examples we take as soft supersymmetry-breaking parameters rfiQ = rfiT = rriB = rrisq 
and At = Ab = A, in the conventions of refs. [^,|^. For comparison, we also show the values 
of the width obtained by replacing^ the amplitude A in eq. @ with its 'improved tree-level' 
expression, A^^'^'^ = Xuhh, and with the effective-potential expression, A^^'^' in eq. (P^j). The 
behaviour of the one-loop-corrected width T{H — > hh) in the (mA,tan/?) plane is illustrated 
in fig. 3, which displays contours of constant width, for the representative parameter choice 
mt = 140 GeV, m,, = 1 TeV, A = /i = 0. 

As a first comment, we observe that there is a region of the (r7iA,tan/5) plane in which the 
decay H ^ hh is kinematically forbidden. At tree level, this region corresponds to | cos 2 f3\ > 
2{m\ + m|)/(5m^m^), which implies mz/'i < rriA < 2m^ and tan/3 > 3. As expected, 
this region is deformed by the inclusion of radiative corrections, in a way which depends on 
mtjfnsq,. ... It is delimited by the thick solid line in fig. 3, and its existence is also evident 
in figs. 2b and 2c. For a given tan/3 in the above range, the forbidden region for H ^ hh 
essentially corresponds to m™" < ruH ^ 2m^°'^, where m™" (m^"^) is the lowest (highest) 
possible value of rriH {nih). We also recall that the small region of rriH ~ -m™", corresponding 
to mA_ ^ 50 GeV, is almost entirely ruled out by the present LEP data |T^. From figs. 2a, 2d 
and 3 we can also see that there is an additional line in the {rriA, tan/3) plane where r{H -^ hh) 
vanishes, due to the vanishing of the amplitude. For the parameter choice of fig. 3, this occurs 
for tan/3 ^ 2 and m^ ~ rriw, corresponding to mn ~ m^*'^. 

The general behaviour of r{H —>■ hh) in the (m^,tan/3) plane is well represented in fig. 3. 
For rriA ^ 2m z, and in the kinematically allowed region, the decay rate depends mildly on 
tan/3, and rapidly decreases with ttia approaching the critical line near mw For m^ ^ 2mz, 
the partial width has a milder dependence on m^ and a stronger dependence on tan/3: the 
largest values are obtained for tan/3 ~ 2—3 and rriA ~ 200—350 GeV. For very large values 
of tan/3 the width becomes negligibly small in comparison with the competing channnels, in 
particular H -^ bb. 

As for the dependence of the corrections on mt, A and fi, in general r{H —>■ hh) rapidly 
increases with increasing mt, A and /i. For example, for m^ = 500 GeV, tan/3 = 1.5 and 
m,sq = 750 GeV, one obtains 

mi = 120GeV,A = /i = 250GeV, ^ F = 0.04 GeV 

mt = 120GeV,A = /i = lTeV, => F = 0.06 GeV 

mi = 180GeV,A = /i = 250GeV, => F = 0.18 GeV 

mi = 180GeV,A = /i = lTeV, => F = 0.37 GeV. 

The comparison between the present one-loop calculations and previous approximations can 
be done by looking at fig. 2. One can see that the 'improved tree-level' result, obtained by using 
the tree-level formula with one-loop-corrected values m,^, m,H and a, can be off by as much as 
a factor of ~ 4. The pure tree- level calculation would be in general in much worse agreement. 
On the other hand, the effective potential result is typically much closer to the full one. In 
particular, the agreement betweeen the two methods is good for tan/3 3> 1 or m,H < 2mt, 2m.t^. 
As expected, radiative correction effects are maximal for tan/3 ~ 1, corresponding to maximal 
top Yukawa coupling for a given top mass. Also, momentum-dependent effects begin to play a 
role only when the top or the stop thresholds are approached. The effect of the top threshold 



^In the evaluation of r(_ff -^ hh), one also needs the masses ?7i^ and uih as functions of the input parameters: 
for those we use the one-loop-corrected expressions in all three cases. To be consistent with this prescription, 
even when evaluating A*^'^'^ we use the one-loop-corrected expression of the mixing angle a. 



is always smooth, but is nevertheless clearly visible in figs. 2a-2d. Besides their effect on the 
proper vertex, the stop thresholds give rise to singularities in Il'fjfj{m'jj), corresponding to a 
breakdown of perturbation theory near threshold: an example is the cusp appearing in fig. 2d, 
whose details should therefore not be trusted. 

To allow a better understanding of the origin of our results, we describe in more detail the 



sources of numerically large corrections. Wave-function renormalization effects in eq. (|T2D are 
in general negligible, also thanks to the use of the one-loop-corrected mixing angle a in the 
definition of the (if, h) basi^. Large corrections come only from the proper vertex Anhh, with 
the imaginary part never very large but sometimes non- negligible, and diagrams of type (II) 
usually give small contributions, unless H is close to a stop threshold. 

We conclude with some comments on the phenomenological implications of our results. To 
consistently examine the effects of these corrections on the H branching ratios, one should also 
compute the remaining partial widths at the one-loop level, which goes beyond the aim of the 
present paper. At present, only approximate computations exist, but some quahtative con- 
siderations are nevertheless possible. For example, one might wonder if the phenomenological 
analyses of H signals at hadron colliders |T0|-|12|, which used the effective potential approxi- 
mation in the computation of r{H -^ hh), are going to be significantly affected. The answer 
is negative: H ^ hh is an important decay mode only for small tan/5 and ttlh ^ 2mt,2mi^, 
and in this region the effective potential approach is a rather accurate approximation. The 
residual small corrections are negligible compared with the uncertainties in the evaluation of 
the production cross-sections. However, one could think that, in the happy event that H is 
discovered, the effects studied in the present paper might be important for a detailed study of 
its properties at a high-energy and high-luminosity e^e~ collider. 
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Appendix 



We give here the exphcit expressions of the coefficients Cjg^g^ and Cyg^g^: 

,2 , 



^Hht2 
^Hb2b2 

'^Hbih 



^ht2i2 



C/itita 



^hhibx 



^hb2b2 



^hbib2 



diim^ cos{a + f3) + 



sin/3 



+ 2stCtBtH 



,t _2 „„„/ , m , "^tsma 



'^22"^z cos(a + /3) + 



sin/3 
(i*2"^| cos(a + /3) + (c^ — sf)BtH 



- 2stCtBtH 



d\im\ cos(a + /3) + 



cos/? 



+ 2sbCbBbH 



3b 2 „„„/„ , ^N , "^fccosa 



'^22"^z cos(a + /3) + 



cos/3 
'^i2"^l cos(q; + 13) + {cl - sl)BbH 



— 2sbCbBbH 



,t 2-1 , ^N mucosa 



diirriz sin(a + /3) — 



sin/3 

,2 



7t 2 .;„^„, , a\ ^tcosa 



(^22^^ sin(a + /3) 



sin/3 



+ 2stCtBth 
- 2stCtBth 



d\^ml sin(a + /3) + (c? - sl)B^ 



th 



d^^m^ sin(a + /?) + 



m^ sin a 
cos/3 



2sbCbB, 



bh 



6^22^1 sin(Q; + /3) + 



cos/3 



+ 2sbCb5| 



bh 



(il2"^| sin (a + /3) - (c^ - 5^)5^,^ 



ml 



^Hhtatt A 



sin 2a 2d\bm\ - Sab^^ 



2 
* _2 , ;r "^b 



sin /3 



%hb.6, = 7 sin 2a 2d„fem^ + 6ab 



cos^/? 



1 / .+ o _ . m? cos^ a ' 



Cfcht„4 = 7 ^afe"^z cos 2a - 6ab 



sin"'/3 



' ^hhbabi, ~ 2 \ "'ab"''Z 



- d^„i,m% cos 2a — 5, 



afe" 



m^ sin^ a ' 
cos^/3 



The conventions for the squark masses and mixing angles, and the symbols (i^^, BqH, Bg^, 
etc., were all defined in [|],0. The symbol 6ab is the Kronecker delta. Notice that the cnhgagb 
coefficients disagree with those reported in ref. [|l^], as already observed in ref. |ll5| . 

The function F^, corresponding to the two-point scalar loop integral, is given by 



'\ 



nil + 1^2 



ms 



1- 



frrii — 7n2 



V m3 



F'^{'mi,m2,m3) = F{mi^m2-,m^)—iT\9{m^—mi—m2) 

The function f(ma,'mb,mc;pi,p2,P3) corresponds to the three-point scalar loop integral 

f{ma,mb,mc;pi,p2,P3) = dx dy pjx^ + ply'^ + {pi - pj - pl)xy 

Jo Jo L 

I / 2 I 2 2\ I / 2 2 I 2 2\ I 2 • T "^ 
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It was studied and computed in ref. [0. For an on-shell decay its explicit expression can be 
written in the form 



f{ma,mb,mc;pi,p2,P3) 



D 



SpiAi^) + SpiAt_) - SpiA^^) - SpiA^_ 



1^2) + (1^3) 



where Sp{x) is the Spence function and 

±pI -ml + ml + Bi 



4± 



B, 



Bi±Ci 



1 



— [PliPl -P2-P3 + 2"\ - "^a - ^b) + (P-i - P2)i^a - ^b)] 



Ci = \Jp\ - 2{ml + niDpl + {ml - m^)^ + ie , 
D = 



ip\ + pi + pi- 2pIpI - 2pIpI - 2pIpI . 

The replacements (1 -^ 2) and (1 — > 3) mean the following: 



1^2 = {pi,p2,P3,ma,mb,mc 
1^3 = {pi,p2,P3,ma,mb,mc 



{p2,P3,Pi,mb,mc,ma) 
{p3,Pi,P2,mc,ma,mb) 



Analogously to the function F'^, the function / develops an imaginary part above thresholds. 
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Figure captions 

Fig.l: The three basic topologies of the diagrams involving top, bottom, stop and sbottom 

exchanges and contributing to the one-loop Hhh proper vertex (g = t, b; a, b, c, = 1, 2). 

Fig. 2: The decay width r{H -^ hh), as a function of rriH, corresponding to the four indicated 
parameter choices. Solid lines correspond to the full diagrammatic calculation, dashed 
lines to the effective potential approach, dash-dotted lines to the 'improved tree-level' 
result. 

Fig. 3: Contours in the (myi,tan/3) plane, corresponding to constant values of T{H -^ hh), for 
the representative parameter choice rrit = 140 GeV, rrisq = 1 TeV, A = fi = 0. The thick 
solid lines correspond to r{H — > hh) = and in particular delimit the region where the 
decay H —^ hh is kinematically disallowed. 
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